Nearly universal crossing point of the specific heat curves of Hubbard models 



OS 

o> 

u 

Of 

< 

o 



D 
i 

CO 

s 

I 

C 

o 
o 



> 

m 
o 

00 
On 

a 

S3 
O 
O 



X 

S3 



N. Chandra,0 M. Kollar, D. Vollhardt 

Theoretische Physik III, Elektronische Korrelationen und Magnetismus, Institut fur Physik, Universitat Augsburg, 

86135 Augsburg, Germany 
(October 28, 1998) 

A nearly universal feature of the specific heat curves C(T, U) vs. T for different U of a general 
class of Hubbard models is observed. That is, the value C+ of the specific heat curves at their 
high-temperature crossing point T+ is almost independent of lattice structure and spatial dimension 
d, with C+/fcs w 0.34. This surprising feature is explained within second order perturbation theory 
in U by identifying two small parameters controlling the value of C+ : the integral over the deviation 



of the density of states N(e) from a constant value, characterized by 5N = Jde \N(e) 
inverse dimension, 1/d. 
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I. INTRODUCTION 

Recently attention was drawn to the fact that in var- 
ious strongly correlated systems the curves of the spe- 
cific heat C(T,X) vs. temperature T cross once or twice 
when plotted for different values of a second thermody- 
namic variable X .El For example, crossing points are ob- 
served for different pressures (X = P) in normalnuid. 
3 He (Ref. 0) and heavy-fermion systems such as CeA^El 
and UBe3D By changing the magnetic field (X = B) the 
same feature is-seen in heavy-fermion-compounds such 
as CeCug-zAlJa and Nd 2 _ :r Ce :c Cu04.tl Crossings of the 
specific heat curves are also observed in the simplest lat=. 
tice model for correlated electrons, the Hubbard model,El 



kcr^ko- 



(i) 



ka 



where ek is the dispersion of a single electronic band, /i 
the chemical potential, and U the local interaction. At 
half filling the curves C(T, U) vs. T always cross at two 
temperatures. This is observed, for example, in the casa 
of the .-model with nearest-neighbor hopping in d = IJjlI 
d = 2 ,EIj .and d — oo |L2l as well as for long-range hopping 
in d = 1 ;E£l for the latter two systems the specific heat is 
shown in Fig. |l|. Furthermore, crossing is found irpi = 1 
when a magnetic field B is changed at constant C/.Ej The 
fact that these crossing points may be very sharp was 
analyzed in Ref. [j], and was traced to the properties of 
certain generalized susceptibilities of the system. 

In the following, we will consider only the crossing of 
the specific heat curves occuring for X = U in the para- 
magnetic phase of the Hubbard model with a symmetric 
half- filled band (n = 1). We will investigate yet another 
observation, namely that for small U the specific heat at 
the high-tcmperature crossing point has practically the 
same value of approximately 0.34fce for all dimensions 
d and dispersions £k, which can be seen also in Fig. [l]. 
This is surprising, because the temperatures at which 
this crossing occurs are very different for different disper- 
sions and dimensions, and because the maximum value 
of C(T, U) and its value at the low-temperature crossing 
point vary strongly, too. It should be noted that the spe- 



cific heat, like the entropy, is a dimensionless quantity 
when expressed in units of ks- 

We denote by T + the temperature at which the curves 
C(T, U) vs. T cross for different values of U. Then the 
specific heat is independent of U at the crossing temper- 
ature T + (U), defined bytl 



dC_ 

dU 



= 



(2) 



T+(U) 



Since we are not concerned with the dependence of T + (U) 
on U, but rather with the crossing point value of C(T, U) 
for different lattice systems, we consider only the limit of 
small U and define 



C+ = lim C(T + (U),U). 



(3) 
(4) 
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FIG. 1. Specific heat for the Hu-bhard model, (a) Exact 
solution for 1/r hopping in d = 1.E3 (bVJterated perturba- 
tion theory for NN hopping in d — tx)£j In (a) T and U 
are in units of the half bandwidth, while for (b) the second 
moment of the (Gaussian) density of states is set to unity. 
At the high-temperature crossing point the specific heat has 
the almost universal value of 0.34fc_B in the limit U — > (see 
arrows) . 
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This weak-coupling crossing point can be calculated 
without approximation using second order perturba- 
tion theory in U. Higher orders in perturbation the- 
ory would be necessary to determine the dependence of 
C(T+(U),U) on U. 

It is the purpose of this paper to show that the value of 
C+ at the high temperature crossing point is almost uni- 
versal, and to analyze the origin of this peculiar feature. 
We calculate C+ for a half-filled band with a symmetric 
density of states (DOS) N(e) = JdkS(e - e k ). We also 
show that the weak dependence of C + on lattice proper- 
ties can be understood by starting from the limit d — oo 
and using expansions in terms of two small parameters. 

This paper is organized as follows. In Sec. |l| we re- 
view the origin of crossing points in the Hubbard model 
and show how to calculate T + and C+ in second order 
perturbation theory, the details of which are contained in 
Appendix [A]. Various non-interacting lattice systems are 
listed in Sec. |ll], and the values of C+ at the high-tem- 
perature crossing point for these systems are presented 
in Sec. 



IV 



Sec. |Vj contains expansions of C+ that re- 
veal the influence of the density of states and the lattice 
dimension. We close with a conclusion in Sec. VI. 



II. 



CROSSING POINTS IN THE SPECIFIC HEAT 
OF THE HUBBARD MODEL 

The entropy per lattice site S(T, U) is given by 



i 

S(T, U) = JdT 



, C (V,U) 



For the Hubbard model S(T, U) approaches a constant 
when T — > oo. Taking the derivative of Eq. (|5|) with 
respect to U we find 



rdT dC(T,U) 

~ J ~T dU ' 

o 



(6) 



Since dC/dU is not identically zero, there must exist 
temperature regions where it has positive and negative 
values. We assume that no phase transitions occur, so 
that dC/dU is a continuous function of U. Then there 
exist temperatures where dC/dU changes sign; at these 
temperatures the curves C(T, U) vs. T cross [see Eq. (||)]. 

There are two such crossing points in the paramag- 
netic phase of the half-filled Hubbard model, as can be 
seen from the sign of dC/dU at very low and very high 
temperatures. For intermediate values of U the specific 
heat of this model shows the following general features: 
C(T, U) starts linear in T at low temperatures and devel- 
ops a two-peak structure, one at temperatures T ~ At 2 /U 
due to spin excitations (t is the hopping amplitude), and 
one at temperatures T ~ U — W due to charge excita- 
tions, where W is the bandwidth. For low temperatures 
the spin exitations become stronger for increasing U, thus 
dC/dU > 0. At high temperatures an increase in U 
pushes out the charge peak and thus increases C(T, U), 
which tends to (a + bU 2 + 0(U 4 ))/T 2 with a, b > 0. 
Hence dC/dU > for both high and low temperatures, 
so that the sum rule [Eq. @] yields an intermediate re- 
gion with OC/dU < 0. There are thus two sign changes 
corresponding to two crossing points. 

To determine the location of the crossing points defined 
by Eqs. (||) and ([|) , we calculate the internal energy per 
lattice site in perturbation theory in U, 

E(T, U) = £ (0) (T) + - U + U 2 E {2) (T) + 0(U 3 ). (7) 



(5) Here (k B =1,0= 1/T) 



E {0 \T) = 2 Jde- 



N(e)e 



exp(/3e) 



(8) 



which is the internal energy for the non-interacting sys- 
tem, {7/4 is the Hartree contribution, and the second- 
order correlation energy is given by (see Appendix ^ for 
details) 



E&)(T) 



80 32 



dx dk dp / dq 



cosh[^a;/3(e k + e p + e p+q + e k + q )] 

cosh(i/3ek) cosh(i/3e p ) cosh(i/3ek+ q ) cosh(i/3e p+q ) ' 



(9) 



where the integrations, e. g. Jdk = fd d k/(2ir) d , run 
over the first Brillouin zone. ,_. 

In the limit of infinite spatial dimensionalj this expres- 
sion can be simplified further. In this case_momentum 
conservation at vertices becomes irrelevant Jl3 so that the 
integrals factorize (see Appendix |A|) 



£ (2) (T) = 



d_0?_ 

'80 32 



dx 



deN(e) 



cosh(ia;/3e) 
cosh(^/3e) 



(10) 



Note that as usual in infinite dimensions the dispersion 



ek enters into one-particle quantities only via the DOS 
N(e). Therefore this expression is much easier to evalu- 
ate numerically than Eq. (||). 

The specific heat C(T, U) = dE/dT has the expansion 

C(T, U) = C (0) (T) + U 2 C {2) (T) + (9(C/ 4 ), (11) 

where 



C^T) = tlde N ^ 



cosh 2 (i/3e) : 



(12) 
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and the function C&(T) can be written as (see Appendix 



C (2) (T) 



32 dfr 



1 

/? 2 [dxJ2[M^0t- 



(13) 



Here the sum runs over lattice sites R m and the functions 
f m (x,/3) are given by 



dk 



exp (i k-R m + ^x/3ek) 
cosh(i/3e k ) 



(14) 



Comparison with Eq. (|T^) shows that in d = oo only 
the local term with R m = 0, i. e. f (x,(3) — Jde N(e) 
cosh(ia;/3e) / cosh(i/3e), contributes to the sum in Eq. 

The crossing point in the specific heat occurs at the 
temperature T + (U) for which C(T, U) is independent of 
U. In view of Eqs. (||), d§|), and (|ll"l), the crossing tem- 
perature T+ in the limit U ~ > is given by the root of 
the equation 

C (2) (T+) = 0, (15) 

and the specific heat at the crossing point, Eq. (jij), is 

C+ = C (0) (T + ). (16) 

These equations will be evaluated for several lattices and 
dimensions that are described in the next section. 



III. MOMENTUM DISPERSION AND DENSITY 
OF STATES 

We consider only one-band systems at half-filling with 
a symmetric density of states on lattices in finite and in- 
finite spatial dimensions. For systems with finite band- 
width we set W/2 = 1, where W is the bandwidth, while 
for infinite bandwidth we use a unit second moment of 
the density of states, i. e. J deN(e) e 2 = 1. 

1. Finite Dimensions. For the linear chain, square 
lattice, and simple cubic lattice, i. e., the hypercubic lat- 
tices in d = 1,2,3, we use the tight-binding dispersion 
£k = — 2t ^2i = i cos hi, \ki\ < 7r, which describes nearest 
neighbor (NN) hopping with amplitude t = l/2d. Fur- 
thermore we study the body-centered cubic (bcc) lattice 
in d = 3 with NN hopping, which can be regarded as a 
subset of the simple cubic lattice with hopping across the 
space-diagonal, so that = — 8t cos(k x ) cos(k y ) cos(fc z ) 
with t = i. For these systems we use perturbation the- 
ory as described in Sec. ||. Finally, for one-dimension- 
al long-rapge hopping t(r) oc 1/r the known interacting 
dispersiontHI can be used instead of perturbation theory. 
The free dispersion is tk = tk, t = with a constant 
density of states. Fig. || shows the various densities of 
states used in d = 1, 2, 3. 




FIG. 2. DOS for several non-interacting systems with 
next-neighbor and 1/r-hopping in dimensions d = 1,2, 3. 

2. Infinite dimensions. We consider first the hyper- 
cubic lattice and generalized honeycomb lattice with NN 
hopping. pFor the hypercubic lattice, the hopping must 
be scaledtil as t = 1/ \2d to obtain a nontrivial limit for 
d — > oo. In this case the density of states becomes a 
Gaussian with unit variance, N{e) — exp (— £i/2)/v / 2~7r, 
whereas for the generalized honeycomb latticeEj the same 
scaling leads to N(e) = \e\ exp (— e 2 ). We also study the 
Bethe lattice with infinite connectivity and semicircular 
density of states, N(e) = f VI - e 2 , |e| < 1. 

Furthermore we can take advantage of the fact that in 
d = oo only the density of states of the non-interacting 
system appears in Eq. (|lfj|). Hence it may be chosen at 
will even when no corresponding dispersion is known. We 
consider three such functions each containing a tunable 
real parameter a > 0. This allows us to study the be- 
havior of C + for a wide range of DOS shapes. The first 
is the "metallic" density of states 



N(e) 



1 



2a 



(i-kD, N<i, 



(17) 



bearing this name due to its finite value at the Fermi en- 
ergy, while the following "semi-metallic" density of states 



N(e) 



1 



e <1, 



(18) 



vanishes at the Fermi energy. Furthermore we employ a 
semi- metallic density of states "with tails" having infinite 
bandwidth and unit variance 



N(e) 



r( c ) 



\e\ a exp(-ce 2 



1 



(19) 



It reduces to the DOS for the generalized honeycomb 
lattice with NN hopping in the case of a = 1. The im- 
portant special case of a constant rectangular density of 
states, N(e) = | for |e| < 1, is contained in Eq. ( |l7| ) for 
a — > oo and in Eq. (18) for a — > 0. Note that in the limit 
a — > oo, the DOS in Eq. (|18|) approaches two (5-peaks, 
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N(e) = + e) + 6(1 - e)). This particular case is of 

interest only because in this case the integrals in Eqs. (||) 
and @ can be calculated analytically. For general a, 
on the other hand, the functions in Eqs. (|l^) and (fig ) 
model typical DOS shapes for finite dimensions. Several 
densities of states used in d = oo are depicted in Fig. [|. 



IV. RESULTS FOR THE SPECIFIC HEAT AT 
THE HIGH-TEMPERATURE CROSSING POINT 



In this section we present results for the specific heat 
C+ at the high-temperature crossing point for the density 
of states discussed in Sec. ED. To calculate C+ accord- 
ing to Eqs. ( 15|) and (|l^) , the integrals appearing in Eqs. 
( |l2| ) and ([lS ) have to be calculated numerically. We de- 
termine them to high precision (typically 10 -8 ) by either' 
Monte-Carlo integration (using the VEGAS algorithmllZl) 
or by a high-temperature expansion (described in Ap- 
pendix |A|) . For several cases both methods were applied 
and yielded the same results within numerical accuracy. 
Typical results for the functions (T) and (T) are 
shown over a wide temperature range in Fig. ^ for the 
linear chain with nearest neighbor hopping and for the 
constant DOS in infinite dimensions. There is a max- 
imum at lower and a minimum at higher temperatures 
in E( 2 '(T), corresponding to the two zeros of C^(T). 
These are the temperatures where the specific heat curves 
cross for [/ — > 0. 

Numerical values for T + and C+ are listed in Table Q. 
In view of the drastically different DOS shapes (see Fig. 
H and [3]) it is quite remarkable that the values of C+ 
for these systems are very similar, ranging from 0.331 
to 0.358. For the DOS with tunable parameter C+ is 
plotted vs. a in Fig. ||. Again, C+ varies only by a small 
amount although the shapes of the DOS change strongly. 
In particular, for the semi-metallic DOS [Eq. Jig)] C+ lies 
between the values for the constant and two 5-peak DOS, 
since these are the limits of Eq. (H8h for a — ► 0, 00. 
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FIG. 3. DOS for several non-interacting systems in infinite 
dimensions. 
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FIG. 4. Second-order contribution in U to the specific heat, 
C^ 2 '(T), for the linear chain (with NN hopping) and for the 
constant DOS in infinite dimensions. The half bandwidth 
is set to 1 in both cases. The inset shows the correlation 
energy/f/ 2 . 



system 


hopping 


d 


c+ 


T+ 


linear chain 


l/r 


1 


0.346994 


0.561816 


linear chain 


NN 


1 


0.355547 


0.705047 


square lattice 


NN 


2 


0.352682 


0.443585 


simple cubic 


NN 


3 


0.348327 


0.358091 


body-cent, cubic 


NN 


3 


0.357578 


0.241221 


hypercubic 


NN 


00 


0.343630 


0.847667 


hyperdiamond 


NN 


00 


0.338411 


0.983569 


Bethe lattice 


NN 


00 


0.340906 


0.480185 


rectangular DOS 


undeterm. 


00 


0.339352 


0.571895 


two 5-peaks DOS 


undeterm. 


00 


0.330857 


1.115358 



TABLE I. Values of the crossing temperature T+ and spe- 
cific heat C+ at the high-temperature crossing point for sev- 
eral d- dimensional Hubbard models. The temperature is given 
in units of the half bandwidth, except for the hypercubic and 
hyperdiamond lattice, for which the variance of the DOS is 
set to unity. 
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FIG 5. C+ for several DOS with a tunable parameter a, 
Eqs. @-(|l9|), plotted vs. a. The upper and lower dotted 
horizontal lines represent C+ for the constant DOS and the 
two 5-peak DOS, respectively. 



4 



These results raise the following questions which will 
be addressed in the next section: (i) Why is C+ at the 
high-temperature crossing point so insensitive against 
changes of the DOS of the non-interacting electrons and 
the spatial dimension? (ii) What determines the (small) 
spread in the values of C+1 (iii) Why is the value of C+ 
at the low-temperature crossing point varying so much 
stronger? 



V. EXPANSIONS OF C 4 



We will see that the influence of the DOS and of the 
crystal lattice on C+ can be understood by expanding 
around the limit of d — oo. In addition, in d = oo the 
dispersion ek enters only via the DOS N(e), so that the 
effect of its form on C+ can be studied by expanding 
in terms of the difference between N(e) and a reference 
DOS N(e). 

I 



A. Influence of the DOS in d — oo 

First we consider infinite-dimensional systems with an 
arbitrary symmetric DOS N(e) with finite bandwidth. 
The DOS is compared to a rectangular-shaped DOS with 
the same bandwidth, N(e) = \ for |e| < 1. Their differ- 
ence is characterized by the quantity 



6N= de\N(e) - N(e)\, 



(20) 



which will serve as an expansion parameter. We will 
expand C+ to lowest order in 5N, making use of the 
known functions fo(x,0), C^(0), C^°'(0) pertaining to 
N(e), as well as the known crossing points for the rect- 
angular DOS, T+ = 1/0+ = 0.13801 and C+ = 0.44046 
at low temperatures, while T+ = 1/0+ = 0.57190 and 
C+ = 0.33935 at the high-temperature crossing point. 
We begin by expanding (0) for small 5N: 



C (2) {0) = C (2) {0) 



2 d 2 
32O0 2 



2 dx[f o (x,0)-f o (x,0)}f o (x,0) 3 



0((5N) 2 ). 



(21) 



Using standard inequalities it can be shown that the neglected terms indeed vanish like (SN) 2 (for fixed < oo). 
Now C^ 2 \0) is expanded around the known crossing point 0+ where C^(0+) = 0. Then the new crossing point 
temperature is determined from the condition 



O = C^{0 + ) = {0+-0+)C^'(0+)+A 
To lowest order the shift in 0+ is hence given by 



2 d 2 ,, 



32 d0 2 



Z dxf o (x,0)f o (x,0Y 



0((SN) 2 



(22) 



/3=/3 4 



50+ =0+-0+= deN(e) A0+(e) + 0((SN) 2 ), 



(23) 



where 



A/3+(e) 



Pi 



\CW(0+) 



B 2 



d0 2 cosh(i^e) 



dx fo(x, 0) 3 cosh(^x0e) 



Finally we expand (0) in — 0+ to first order and 
evaluate it for — 0+ + 50+. Thus the specific heat at 
the crossing point C+ is obtained to lowest order in 5N 
as 



C+ = C+ + JdeN(e)AC+(e) + 0{(5N) 2 ), (25) 

-l 



where 



(24) 



P=0+ 



AC+(e) 



2 cosh^(i/3 + e) 
+ C^'(0+)A0+(e) 



CW(0+). (26) 



The last two equations show how the first order effect on 
C+ of a deviation of N(e) from a rectangular shape can 
be determined by a single integration. 

The numerical evaluations of the function AC+ (e) are 
plotted in Fig. ^| for both crossing points. The ampli- 
tudes of the function AC+(e) corresponding to the low- 
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and high-temperature crossing points are seen to differ 
greatly, i. e. by a factor of about 40. This implies a 
much greater sensitivity towards changes of the DOS and 
the dimension of the value of C+ at the Zow-temperature 
crossing point. 

It is also clear that C+ is not entirely universal at the 
high-temperature crossing point, since for a general den- 
sity of states N(e) the integral in Eq. (|25|) does not van- 
ish. We can in fact estimate the maximum value of the 
shift <5C+ = C+ — C+ for arbitrary N(e) with finite band- 
width, using simple integral inequalities: 



Hence we have shown that for a small change SN in the 
density of states the variation of C+ at the low-tempera- 
ture crossing point is large, while at the high-temperature 
crossing point it is small and well-described by the first 
order correction in SN. This gives quantitative answers 



to the questions posed at the end of Sec. IV 



B. Lattice effects in dimension d = 1, 2, 3 



\SC+\ < min ^oi, a 2 ■ max [iV(e)] 



0((<«V) 2 ), (27) 



with ai = max < e <i |AC+(e)|, a 2 = f_-.de |AC+(e)|. At 
the high-temperature crossing point we find a\ — 0.02268 
and a 2 = 0.02200, i. e. to order 0(5N) we have \6C+ \ < 
C+, which is the reason for the insensitivity of C+ to 
changes in N(e). Furthermore the predicted range of val- 
ues 0.339 ± 0.023 corresponds well to the observed range 
of 0.331-0.358; see Table @. On the other hand, at the 
low-temperature crossing point we find ai — 0.9330 and 
a 2 = 0.7727, so that \SC + \ w C + . Hence C + is indeed 
not confined to a small interval in this case. 

To check the validity of this expansion we applied Eq. 
( p5j ) to several infinite-dimensional systems at the high- 
temperature crossing point. Results for the Bethe lat- 
tice, as well as for the metallic DOS [Eq. @] and the 
semi-metallic DOS, [Eq. (fl8|)j_a,t several values of the pa- 
rameter a are given in Table _Q As expected we find that 
the lowest-order approximation in Eq. ( |25| ) describes the 
behavior of C + very well if the deviation from the rectan- 
gular DOS is not too large. The difference between the 
exact value of C+ and the approximate value C+ + SC+ 
is due to corrections of order (SN) 2 and hence is typi- 
cally an order of magnitude smaller than the first order 
correction SC+. 
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FIG. 6. Weight function AC+(e) determining the shift in 
C+ at the low- and high-temperature crossing points for a 
Hubbard model in d — oo with finite bandwidth, according 
to Eq. (bq). The half bandwidth is set to unity. 



To determine the influence of the lattice dimension on 
the value of C+ we study JJae hypercubic lattice with NN 
hopping t and the scalingE-3 t = 1 /y2d and expand C + in 
1/d. Since the large variation of C+ at the low-temper- 
ature crossing point can already be understood from the 
results of the previous subsection, we will perform this 
calculation only for the high-temperature crossing point. 

As discussed in Sec. |l[ for d = oo only the local 
term involving fa(x,j3) remains in the lattice sum in 
Eq. (|l3|). The 1/d-corrections are given by the contri- 
bution from the 2d nearest neighbors vectors such as 

R 1 = (i,o,o,...),y 

fm{x,(3Y = / (x,/3) 4 + 2df 1 (x,(3) i + ■■■ 



= ./o(.t,/3) 4 + 



2d 



2d_ 



fo(x,/3) 



0{d- 



(28) 



where a partial derivative w.r.t. x was employed in the 
second line to remove the factor cos(k-Ri) from fi(x, j3). 
Hence the corrections to C+ in order l/d can be calcu- 
lated entirely from fo(x,(3). This function contains an 
integral over the P-OI3 N(e) only, which in turn must be 
expanded in 1/d.Ej Then the finite-ci correction to T + 
and C+ can be calculated by expanding Eqs. ( |l5j ) and 
( |l6| ) to first order in l/d, with the following result for 
the high-temperature crossing point: 



C+ = 0.343630 + 0.013599 - + 0(d" 2 ) 



T, 



0.847667 + 0.082650 - + 0{d~ 2 ) 



2d 



(29) 
(30) 



Here the temperature scale has been reset to our previous 
choice of unit half bandwidth in finite dimensions (i. e., 
the extra factor l/\/2d must be omitted in order to re- 
cover T_|_ in d = oo). Again this expansion compares very 
well with the numerical results; see Table 111. Note that 
the coefficient of d^ 1 in Eq. ( |29| ) is already much smaller 
than C+ in d = oo, which is the reason for the insen- 
sitivity of C+ to lattice effects, such as Brillouin zone 
shape and momentum conservation. We expect that for 
longer-range hopping the smallness of the deviations in 
the value of C + may equally be traced to 1/ci-corrections. 
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VI. CONCLUSION 



why C+ has an almost universal value at the high-tem- 
perature crossing point. 



For many correlated electronic systems the curves of 
the specific heat vs. temperature obtained for differ- 
ent values of a second thermodynamic parameter X are 
known to intersect n For the Hubbard model at half- filling 
the specific heat curves for different values of the Hub- 
bard interaction U cross twice, the crossing point at high- 
temperatures being remarkably sharp up to intermediate 



values of U. We observed that C + , the value of the spe- E(T, U) = T 
cific heat at this crossing point in the weak-coupling limit, 
is practically the same for a several different lattice sys- 
tems. We analyzed the origin of this conspicuous feature 
by calculating C+ in perturbation theory in U. We found 
the values of C+ at the high-temperature crossing point 
to occur in very small interval, i. e. C+ f=a 0.34 is indeed 
almost independent of dimensionality, crystal lattice, and 
energy dispersion. This is not the case for the crossing 
point at low temperatures, where C+ varies on a much 
larger scale. 

Qualitatively, the reason for this difference can be 
traced to the relevant energy scales on which C(T, U) 
varies. At high temperatures, the energy scale for T 
is essentially determined by the bandwidth, i. e. by the 
hopping amplitude t in the dispersion et- At low temper- 
atures, on the other hand, the generation of low-energy 
excitations (which are responsible for the strong enhance- 
ment of the low-temperature specific heat and dC /dU > 
0) leads to a renormalized energy scale t — ► t c g <C t. 
The first maximum in C(T, U) occurs at a temperature 
that is of the order of i e ff (see Fig. ||). As a consequence, 
the first sign change in dC/dU is also linked to i e ff, so 
that the intersection of C(T, U) and C(T, 0) at low tem- 
peratures does not occur at any predetermined value. In 
contrast, the second sign change in dC/dU is determined 
only by energy scales that also appear in the non-inter- 
acting system, leading to a nearly universal value for the 
high-temperature intersection of C(T, U) and C(T, 0). 

To gain a more quantitative understanding we iden- 
tified two small parameters which determine the cross- 
ing point values C+. The starting point for expansions 
in these small parameters is the limit of infinite dimen- 
sions (d = oo). (i) For d — oo the dependence of C+ on 
the shape of a DOS N(e) with finite bandwidth is well 
described by the first-order correction in the parameter 
SN = fde\N(e) — i|. This parameter is a measure of 
the difference between N(e) and a constant rectangular 
DOS. It turns out that at the high-temperature crossing 
point this correction is small for almost all DOS, while 
it is large at the low-temperature crossing point, (ii) For 
hypercubic lattices in dimensions I < d < oo the value 
of C + may be obtained by an expansion around d = oo 
in powers of 1 jd. At the high-temperature crossing point 
the value of C+ of the e?-dimensional system is already 
accurately determined by the first-order correction in 1 jd 
even for d as low as d = 1, due to the smallness of the 
prefactor of this term. These expansions show in detail 



APPENDIX A: INTERNAL ENERGY AT WEAK 
COUPLING 

The internal energy per lattice site E is given byEl 

£k + i£ CT (k,iw„) 



dk e 



G (iwn , k) 1 - E CT (k, iw r , 



(Al) 

where w„ denotes fermionic Matsubara frequencies, r\ — > 
+ , and h = 1. The non-interacting Green function is 
Go(iuj n , k) _1 = iuj n — (ek — jtx), and S (T (k, ito n ) is the self- 
energy for spin a. In the paramagnetic phase the sum 
over spins just gives a factor of two, and the spin indices 
on the self-energy can be dropped. 



DOS 


a 


6N 


C+ ,approx 


C+, exact 


diff. 


rectangular 







0.339352 


0.339352 





Bethe lattice 




0.231 


0.341391 


0.340906 


0.14% 


metallic 


4 


0.267 


0.341271 


0.340444 


0.24% 


metallic 


6 


0.207 


0.340482 


0.339966 


0.15% 


metallic 


8 


0.169 


0.340091 


0.339743 


0.10% 


metallic 


10 


0.143 


0.339871 


0.339621 


0.07% 


metallic 


12 


0.124 


0.339735 


0.339552 


0.05% 


semi-metallic 


0.1 


0.070 


0.338349 


0.338381 


0.02% 


semi-metallic 


0.5 


0.296 


0.335393 


0.335880 


0.14% 


semi-metallic 


1 


0.500 


0.333279 


0.334253 


0.29% 


semi-metallic 


2 


0.770 


0.331629 


0.332752 


0.35% 



TABLE II. Comparison of the approximate results for C+ 
at the high-temperature crossing point obtained from the ex- 
pansion in SN [Eq. (^)] with the exact values for several 
DOS in d = oo. The expansion is controlled by the parameter 
SN which measures the "distance" of N(e) to the rectangular 



DOS, see Eq. (|20 
the DOS of Eqs. (1 



'Metallic" and "semi-metallic" refer to 
and (|li|), respectively. The last column 
shows the difference between the approximate and the exact 
results in percent. 



lattice 


d 


C+,approx 


C+, exact 


diff. 


hypercubic 


oo 


0.343630 


0.343630 





simple cubic 


3 


0.348164 


0.348327 


0.05% 


square lattice 


2 


0.350430 


0.352682 


0.64% 


linear chain 


1 


0.357229 


0.355547 


0.47% 



TABLE III. Comparison of the approximate results for C+ 
at the high-temperature crossing point obtained from the ex- 
pansion in 1/d [Eq. (p9|)] with the exact values for hypercubic 
lattices in d dimensions. The last column shows the difference 
between the approximate and the exact results in percent. 
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For a symmetric DOS the chemical potential fi at the new chemical potential fixed at 0. Up to second 
half filling is given by U/2 for all temperatures due to order in U S is given by only one Feynman diagram 
particle- hole symmetry. It is useful to define shifted where the lines now represent Go- We write £(k, iuj n ) = 
functions G^ 1 = G^ 1 - U/2 and t = £ - U/2, with U 2 ■ a(k,ioj n ) + 0{U 3 ), with 
I 



C5=-T 2 ^ Jdk 2 Jdk 2 Jdk 3 



x G (k 2 ,iw/)G (k 3 ,icj ; + iv m ) G (k 4 ,iw„ + w m ) 53 - k 2 + k 3 - k 4 - K). 



(2) 



K 



where v m denote bosonic Matsubara frequencies and the sum is over reciprocal lattice vectors K. Expansion of the 
internal energy [Eq. flAlD ] in powers of U yields Eqs. (0), (^), and 



E {2) {T)=T dkJ2 



lUJ n - e k 



a(k,iuj n ) + 2e k 



<r(k, iuj n 



iu n - e k 



(3) 



The frequency summations are carried out as usual Jla with the result 
<r(ki, iu> n 



iu) n - e kl 



= - dk 2 dk 3 dk 



smh[^/3(e kl - e k2 + e k3 - e k J] ^ 

77TZ — n ■ 1^ 5 ( kl - k 2 + k 3 - k 4 - K), 

£ki - £k 2 + £k 3 - e kJ ll i= i cosh(^Pe ki ) K 



(4) 



Here we eliminate the energy denominator using the identity sinh(i/?y) = i/3 J Q dx cosh(ix/3y). Next we rewrite 
the second term in square brackets in Eq. (||) as a formal derivative 2e kl 9/9e kl of Eq. We make the integrand 
symmetric in all e ki by shifting k 2 and k 4 by Q, where e k+ Q = — e k . Then the derivative can be replaced by 2f5d/d[3 
because only products /3e ki appear. Taking (3 inside the derivative we finally arrive at 



£ (2) (T) = -|^ Jdx Jdk! Jdk 2 Jdk 3 Jdki 

o 

cosh[±x/3(e kl + e k2 + e k;3 + e k4 



cosh( i/3e kl ) cosh( i/3e k2 ) cosh( i/3e k3 ) cosh( i/3e k4 



k 4 - K) 



(5) 



K 



In d = 1,2,3 this simplifies to Eq. ([}]), while for d = oo the ^-function can be omittedO and the numerator can be 
replaced by nt=i cosh(^xPei), leading to Eq. 

The integrals in Eq. fla) also factorize if we express momentum conservation as a sum over lattice vectors R m , 



53 <K k i - k 2 + k 3 - k 4 - K) = 53 ex P M k i - k 2 + k 3 - k 4 - K) • R m ] 



(6) 



K 



leading to Eqs. ( |13| ) and (J14|). Since the functions f m {x,0) in Eq. ( |14[ ) cannot be calculated in closed form for general 
e k we employ a high-temperature expansion, which yields 



f m (x, (3) = 53 ^ £7„(*±*) /„(R m ), 
I„(R m ) = y dk (e k )" exp (i k • R m ), 



(7) 
(8) 



where E n {x) are Euler polynomials. The remaining Brillouin zone integrals are calculated as follows. For the d-di- 
mensional hypercubic lattice the lattice vectors are R m = (mi, . . . ma) with integer m„ so that for NN hopping t we 
have 



/£ c (R m ) = (-*)" 53 



??■ 

m, . . .rid 



n 



7T 

/" ^ (2 cos fci) n< cos(fc l |m l |) 



(9) 
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The integral in square brackets equals if < rii — \mi \ = 2r^ with integer r^, and zero otherwise. Hence 



I^(R m ) = (-t) n n\ J2' II 



r!,...,r d i- 



L n\(ri + \m,i\)V 



(10) 



where the sum is restricted to 2r 1 + . . . + 2r^ = n — M hc , with M hc = Yli=i \ m i\- Note that Z^ c vanishes if M hc > n 
or if n + M hc is odd. For the bcc lattice in d = 3 we use the hypercubic lattice basis R m and the dispersion given in 
Sec. 111. We obtain 



j£ cc (R m )=H)»n 



7T 

/dk 
(2cosfci)™ cos(ki\mi\) 



(-t) n U 



(11) 



with Ti defined as above. In particular I^ cc vanishes if Af bcc = max; \nii\ > n or if n + \rrii\ is odd. 

The advantage of the present high temperature expansion is that terms with M > n vanish when expanding to 
order (3 n , so that the lattice sum terminates at m.; = n. In addition, lattice symmetries can be used to further reduce 
the computational effort. This makes the calculation feasible even in dimension d — 3 where numerical Brillouin zone 
integration is usually difficult. At the high-temperature crossing point 0(/3 80 ) is typically sufficient to obtain T + and 
C+ to an accuracy of 10~ 8 . 



* Present address: Sektion Physik, Ludwig-Maximilians-Universitat Miinchen, Theresienstr. 37, 80333 Munchen, Germany. 
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